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Abstract 

The fluid dynamics video presented here outlines recent advances 
in the simulation of multiphase cellular blood flow through the direct 
numerical simulations of deformable red blood cells (RBCs) demon- 
strated through several numerical experiments. Videos show RBC de- 
formations in variety of numerical simulations, relative viscosity of a 
suspension of RBCs in shear, and the cell-depleted wall layer for blood 
Hagen-Poiseuille flow. 

The animations presented in the submission are results from the direct 
numerical simulation of deformable RBC suspensions using a hybrid lattice- 
Boltzmann spectrin-link (LB-SL) method[lj that is coupled with a particle- 
particle and particle-boundary contact model of MacMeccan et al.[2]. The 
implementation is based on the parallel framework of Clausen et al.(3] has 
also been optimized to run on as many as 65,536 cores of the IBM Blue 
Gene/P architecture. In order to simulate high volume fraction (hema- 
tocrit) simulations of deformable particles in wall-bounded or cylindrical 
domains, a seeding method is used. In this procedure, the particles are 
initialized at 30% of their final size and grown to their full size rigidly, al- 
lowing particle-particle and particle-wall interactions using the lubrication 
and contact model previously referenced. Each RBC is composed of a trian- 
gulated spectrin network consisting of iV=613 nodes. Physically, RBCs are 
biconcave in shape with a maximum diameter of 8 urn enclosed by an elastic 
membrane with an effective elastic shear modulus of 5.7 x 10 -3 dyn cm -1 [4J 
and a bending stiffness of 2.2 x 10~ 12 dyn cm |5], which surrounds a liquid 



hemoglobin with a viscosity of 6 cP. The RBCs are suspended in blood 
plasma with a viscosity of 1.2 cP. The relative non-dimensional elastic pa- 
rameter is the elastic capillary number, Ca<^ = H"ya/G m where [i is the fluid 
viscosity, 7 is the local shear rate, a a length scale based on the particle size, 
and G m is the elastic shear modulus. 

Isolated RBC simulations are first presented. These simulations include 
a numerical experiment chosen to mimic an optical tweezer experiment. 
For isolated RBCs in shear, two different dynamic regimes are then shown 
demonstrating the tumbling versus tank-treading phenomena. A parachut- 
ing RBC is then animated in a microvessel sized rigid tube. The flow in 
that simulation is driven by an axial body force used to imitate a constant 
pressure- gradient . 

Dense suspensions of RBCs follow starting with a demonstration of the 
non-Newtonian effect known as shear thinning is then given by results of 
simulations using the Lees-Edwards boundary condition (LEbc) for a sus- 
pension of spectrin-link based RBCs. As the hematocrit (or volume frac- 
tion) of RBCs increases, so does the relative viscosity. This can be shown 
using the same simulations used to demonstrate shear thinning with the 
LEbc. The cell depleted wall layer is then highlighted through an animation 
of a microvessel sized tube with a focused frame near the tube wall. En- 
semble averaged hematocrit and velocity are also given to demonstrate the 
cell depleted layer thickness and the blunting of the velocity profile due to 
the presence of RBCs. The accumulation of shear stress on the surface of 
platelets is one way to model "blood damage" through blood damage index 
(BDI) models. By simulating a dense suspension of RBCs with a realistic 
number of platelets we can monitor the effects of RBC-platelet interactions 
and the effect on the platelet shear stress. 
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